Infrared catastrophe and tunneling into strongly correlated electron systems: Exact 
solution of the x-ray edge limit for the ID electron gas and 2D Hall fluid 
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In previous work we have proposed that the non-Fermi-liquid spectral properties in a variety of 
low-dimensional and strongly correlated electron systems are caused by the infrared catastrophe, 
and we used an exact functional integral representation for the interacting Green's function to map 
the tunneling problem onto the x-ray edge problem, plus corrections. The corrections are caused by 
the recoil of the tunneling particle, and, in systems where the method is applicable, are not expected 
to change the qualitative form of the tunneling density of states (DOS). Qualitatively correct results 
were obtained for the DOS of the ID electron gas and 2D Hall fluid when the corrections to the x-ray 
edge limit were neglected and when the corresponding Nozieres-De Dominicis integral equations were 
solved by resummation of a divergent perturbation series. Here we reexamine the x-ray edge limit 
for these two models by solving these integral equations exactly, finding the expected modifications 
of the DOS exponent in the ID case but finding no changes in the DOS of the 2D Hall fluid with 
short-range interaction. We also provide, for the first time, an exact solution of the Nozieres-De 
Dominicis equation for the 2D electron gas in the lowest Landau level. 

PACS numbers: 71.10.Pm, 71.27.+a, 73.43.Jn 



I. INTRODUCTION 

In a previous paper^i we proposed a connection 
between anomalies in the tunneling density of states 
(DOS) at the Fermi energy of a wide variety of low- 
dimensional and strongly correlated conductors, and 
the infrared catastrophe. The latter is a well-known 
singular screening response of an ordinary metal to the 
sudden appearance of a localized potential, in this case 
produced by an electron added to the system during a 
tunneling event. Systems where we expect this connec- 
tion to apply include all ID electron systems 
the 2D diffusive metal 9 ! 10 : 11 ! 12 . 13 ! 14 : 15 ! 16 ! 17 ! 18 ! 19 ! 20 
and Hall fluid, 21 ! 22 . 23 ! 24 ! 25 ! 26 . 27 ! 28 : 29 ! 30 

and the edge of the confined Hall 

flujj 31.32.33.34.35.36.37.38.39.40.41.42.43.44.45.46.47.48.49.50.51.52.53 

We argued that in such systems, the accommodation 
of a new electron added during a tunneling event is 
frustrated by the low dimensionality, or the localizing 
effects of a magnetic field or disorder, or both. In these 
cases the tunneling problem is similar to the x-ray edge 
problem. 

A mapping between the two is made via an exact 
scalar functional integral representation for the interact- 
ing propagator, which replaces it by a Gaussian average 
of noninteracting propagators for electrons in the pres- 
ence of potentials cf>(r, r). We then singled out a danger- 
ous field configuration 



Mr,r) =U(t)Q(t -t)Q(t), 



(1) 



which would be the potential produced by an electron 
added to the origin at time r = and removed at a later 
time To , if it had an infinite mass. Here U (r) is bare 
electron-electron interaction. In Ref. |l| this special field 
configuration was treated by resumming a divergent se- 
ries, and fluctuations about </> xr (r, t) were ignored, yet 



qualitatively correct expressions for the DOS were ob- 
tained. 

To obtain quantitatively correct results it will be nec- 
essary to go beyond this "perturbative" x-ray edge limit. 
In Ref. |64| we proposed and investigated a functional 
cumulant expansion method that includes field fluctua- 
tions away from xr (r, r), and treats field configurations 
close to xr (r,r) perturbatively as in Ref. Q. Although 
the improved method yields the exact DOS exponent for 
the important Tomonaga-Luttinger model, calculable by 
bosonization, we do not expect it to be generally exact 
in ID. (Furthermore, the method fails in the presence of 
a strong magnetic field because of the ground state de- 
generacy.) In this paper we neglect fluctuations about 
^xr(r, t) but treat that field configuration exactly (in the 
54r6iaWiTt.5taus6ffteiasar«flptotic limit). This is accomplished 
by finding the exact low-energy solution of the Dyson 
equation for noninteracting electrons in the presence of 
<Axr(r, t), which we refer to as the Nozieres-De Dominicis 
equation. We carry out this analysis for the ID electron 
and 2D Hall fluid, both with short range interaction. To 
this end we obtain, for the first time, an exact solution of 
the Nozieres-De Dominicis equation for the 2D electron 
gas in the lowest Landau level. 



II. FORMALISM 

We calculate the tunneling DOS by analytic continua- 
tion of the Euclidean propagator 

G(r f( r f ,r i( Ti,To) = -<T^ f (r £> 7tt)^(ri,0)) H , (2) 

where H = Hq + V is the grand-canonical Hamiltonian 
for the D-dimensional interacting electron system, with 



E 



n 2 



Mr) 



2 



and 



III. X-RAY GREEN'S FUNCTION 



V=% I d D rd D r' 5n(r)U{r ~r')6n{r'). 



Here II = p + |A, with A the vector potential (if a 
magnetic field is present), and 



Sn(r) 



^^(r)V> ff (r) -n (r) 



(3) 



is the density fluctuation operator. Ho is the Hamilto- 
nian in the Hartree approximation, and v(r) includes any 
single-particle potential along with the Hartree interac- 
tion with the self-consistent density no(r). 

After an exact Hubbard-Stratonovich transformation, 
the interacting Green's function can be written as func- 
tional integral over a scalar field 0, 

x /^[ ( /»]e" 4 ^ t/ "^.g(r f a f ,r i a i ,r o |z0 xr + 0), (4) 



where N= (Texp(— L dTV))^ 1 is a constant, indepen- 
dent of ro, 



(•5) 



e-i/^-V 

is a measure normalized according to J Dfj,[<j>] — 1, and 
g{r i (7 l ,v i c7 i ,T {) \(j)) 

is a noninteracting functional of </>. 

The dangerous field configuration xr (r,-r), which it- 
self depends on the parameters r;, rf , and r appearing 
in @, has been given in Ref. Q]. For the case the tun- 
neling DOS at point ro of interest here we have definition 
Q , which is the potential that would be produced by the 
added particle in J5J if it had an infinite mass. Fluctua- 
tions about ^> xr account for the recoil of the finite-mass 
tunneling electron. 

In the x-ray edge limit, we ignore fluctuations about 
</> xr , in which case 

G(rfcr f ,rif7i,ro) « Afg(r{<Tf , r-^u rol^xr)- (7) 

Eq. Q defines the interacting propagator in the x-ray 
edge limit. The local tunneling DOS at position ro is 
obtained by setting rj = rf = ro and o\ — Of — <7q, and 
summing over <7q . In the remainder of this paper we will 
evaluate for the ID electron gas and the 2D Hall fluid, 
with a short-range interaction of the form 



U(r) = \6(r). 



(8) 



The quantity g(rf<7f , ri<7i, ro|«0 xr ) required in (J7J) is re- 
lated to the Euclidean propagator, 

G xr (rar, r a r ) = {Te -J ^Sn )o < 

according to 

p(rf<Tf,riO-i,To|#xr) = G xr (rfa f T ,ri(TiO) Z^(t ), (9) 
with 

Z xr (r ) = (Te--fo dT I dDr ^ r ^ 5n ( r ^) . (10) 

We refer to G xr (rcrr, r'a'r') as the x-ray Green's function 
which describes noninteracting electrons in the presence 
of a real- valued potential xr (r, t). It satisfies the Dyson 
equation 

G xr (rcrr, r'or') = G (rcr, r'cr, r - r') 

+ J d D fdf G (r<r, fa,T-f) </> xr (f, f) G xr (f erf, rVr'). 

(11) 

Here we have used that fact that the x-ray Green's func- 
tion is diagonal in spin. For a calculation of the DOS we 
use the form Q), in which case l|llf) becomes 

G^{yot,t' or 1 ) = G (r<r,rV,T - r') 

/•TO 

+ X dtG (va,r a,T-t)G w {r at,r'aT'), (12) 
Jo 

where we have assumed the short-range interaction (JSJ). 

By using the linked cluster expansion and coupling- 
constant integration, Z xr can be shown to be related to 
the x-ray Green's function by§£ 

Z xr (T ) = e noXTo e- X ^^o^Io°^ Gi(iWT,ro<rr+) ! ^ 

where G xr (r<7T, r' 'err') is the solution of l|12|) with scaled 
coupling constant £A. 

There is no ro dependence in the DOS for the trans- 
lationally invariant models considered here and we can 
take r = 0. 



IV. ID ELECTRON GAS 

G xr (0aT, Oar') was calculated exactly in the large ro, 
asymptotic limit for the 3D electron gas in zero field 
by Nozieres and De Dominicis^ Their result is actually 
valid for arbitrary spatial dimension D if the appropriate 
noninteracting DOS is used. 

We take the asymptotic form of the noninteracting 
propagator as 

Go(r)w-P— , with N = — . (14) 

T TTVF 
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FIG. 1: The 12 possible time orderings k = 1,2, ••■,12 of 
G X r(0r, Or'), ro is assumed to be nonnegative. 



iVo is the noninteracting DOS per spin component at ep, 
and P denotes the principal part. The solution of 1|12|) 
for this model with r = r' = is 



G xr (r ) = -N cos(6\) 



P y —L + nsm(5 x )6(T ) 



and 



ro 



^xr(To) = — 
\T0 



2(W) 



(15) 



(16) 



the corresponding integral equation for the Hall fluid. We 
were able to guess the exact analytic solution, aided by 
perturbation theory and by numerical studies carried out 
by expansion in a plane- wave basis followed by matrix in- 
version. 

We assume the system to be spin-polarized and spin 
labels are suppressed. In the Landau gauge A = Bxe y , 
the noninteracting propagator in the |r| 3> lu^ 1 limit is 



G (r,r',r)= T(r,r') \v - 6(r)], 



(20) 



where v is the filling factor satisfying < v < 1, and 



r( r ,r') = l p c -lr-r'| 2 /4i 2 c -j(x+x')(y-y')/2e 2 _ (21) 



First consider the case where r; = rf = ro. We can let 
r = without loss of generality and at the origin H12|) 
reduces to 



G xr (0r,0r') 



v - 6(t - t') 



+ 7 / dt[v-@{T-t)] G xr (0i, Or'), (22) 
Jo 



where 



where 5\ is the scattering phase shift of the electrons 
caused by the potential </> xr given by 



5\ = arctan(iVo7rA) 



(17) 



and a is a short time cut-off on the order of the Fermi 
energy. 
Thus 



G(t ) « 3(T |# xr ) 



to 



l+2<5 A /7T+2((5 A /7r) 2 



(18) 



2ir£ 2 



(23) 



is an interaction strength with dimensions of energy. 

The time arguments of G xr (0r, Or') on the left side 
of Eq. (|2"2")l can assume the 12 possible orderings k = 
1,2,- •• ,12 defined in Fig. ^ the right side produces 
terms with two or more different orderings k' , k" , ■ • • . We 
therefore seek a solution of the form 



which gives a DOS in the x-ray edge limit as 



N(e) 



,25 A /7r+2(<5 A /T r ) 2 



(19) 



By expanding the exponent in (|19|l in powers of the cou- 
pling paramenter A one recovers the perturbative x-ray 
result of Ref. 0. 



V. 2D HALL FLUID 

Unlike the low-energy Dyson equation for the ID elec- 
tron gas, which is solvable by Hilbert transform tech- 
niques, there are no standard methods available to solve 



G xr (0r,0r') =^A fc ^ fe (r,r')/ fc (r) 



(24) 



where Wk(r, t') is unity if r and r' have ordering k and 
zero otherwise; an explicit form for Wfe(r, r') is given 
in Appendix [S] The functions fk (r) are chosen to re- 
flect the fact that an electron accumulates an additional 
phase 7 Ar while in the presence of (/> xr for a time Ar, 
whereas a hole acquires a phase —7 Ar. The 12 unknown 
coefficients Ak (which depend parametrically on ro and 
r') are obtained by solving the 12 linearly independent 
equations resulting from the decomposition of (|22|l into 
distinct time orderings k = 1, 2, • • • , 12. The result is 



4 



G xr (0r,0r') = 



- 1) (Wi + W 5 ) + ve-^ ro (W 2 + W 6 ) + {v - l)^ 7 ""^ W 3 



2tt£ 2 \l - v + ise-f T ° 

+ vWn+ve-~< T Wi2 . (25) 



As a side note, the solution for general r; and rf is obtained using the same method and when both r and r' are in 
the interval (0, tq), the result is 



G„(rr,rV) = [u-Q(t - r')} [r(r, r') - 2n£ 2 T(r, 0) T(0, r')] 



+ 27r^ 2 r(r,0)r(0,r') 



the other cases follow similarly. At the origin (|26|l reduces to 



(v - 1) 6(r - r') + ise-'< T <> 9(r' - r) 



1-U+ ve~~< T ° 



-7(t-t'). 



G xr (0r, Or') 



2tt£ 2 



(p - 1) 6(r - t') + i/e- 7To 6(r' - r) 
1 - ^ + ^e-T" 1 "" 



-7(r-r') 



(26) 



(27) 



Finally 

^(0ro,00) = ^( i _; + -^ 7ro )e-^. (28) 

Using Eq. (|T3J) we obtain 

Z XI = e' y7To (l - i/ + i/e^ ), (29) 

therefore 

g(rQa ,r a ,T \i<f> KI ) ■ ( 30 ) 

This is identical to what we obtained in Ref. pj. The 
tunneling DOS is therefore 

N(e) = const x S (e - [l-u]<y). (31) 



VI. DISCUSSION 

In this paper, we have carried out an exact treatment 
of the x-ray edge limit introduced in Ref. 0, for the 
same models considered there. Whereas the ID elec- 
tron gas result (jT§|l would be expected, the DOS of the 
2D Hall fluid remains gapped as in Ref. Q]. A gener- 
alization of our method that accounts for fluctuations 
about xr , and that can be used in a magnetic field, will 
be needed to recover the actual pseudogap of the Hall 

fl mc J 21.22.23.24.25.26.27.28.29.30 



APPENDIX A: TIME ORDERING FUNCTIONS 



Let 



W x {r, 


T ) 


= Q(-t) e(-r') 6(r - t') 


W 2 (t, 


T ) 


= 9(-r) e(-r') 6(r' - r) 




T ) 


= W(t) W(t') 9(t - t') 


W 4 (r, 


T ) 


= W(r) W(t') 6(t' - r) 


W 5 (t, 


r ) 


= e(r-ro)0(T'-T )0(r 


W e (T, 


r ) 


= e(r-ro)e(T'-T )0(r 


W 7 (t, 


r ) 


= 1T(t)0(-t') 


W 8 (t, 


r ) 


= 0(r - r ) O(-r') 


W q (t, 


T ) 


= 0(-t)TT(t') 


W w (t, 


T ) 


= 9(r - tq) W(t') 


W n (T, 


T ) 


= O(-t)O(t'-t ) 


W 12 {t, 


T ) 


= 1T(t)O(t'-t ), 



where ©(<) is the Heaviside step function and W (with 
no subscripts) is the a window function, defined as 



W = Q(t -t)Q{t). (Al) 



■5 
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